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Abstract 


taring the passage of a cold front in March 2002. bio-optica! properties examined in 


coastal waters impacted by the Mississippi River indicated westward adveetive flows and 
increasing river discharge containing a larger nonalgal particle content contributed 
significantly to surface optical variability. A comparison of seasonal data from three 
cruises indicated spectral models of absorption and scattering to be generally consistent 
with other coastal environments, while their parameterization in terms of chlorophyll a 
concentration (Chi) showed seasonal variability. The exponential slope of the colored 
dissolved organic matter (CDOM) averaged 0.0161 ± 0.00054 nm’ 1 , and for nonalgal 
absorption it averaged 0.01 1 nm' 1 with deviations from general trends observed due to 
anomalous water properties. Although the phytoplankton specific absorption coefficients 
varied over a wide range (0.02 to 0.1 m 2 (mg Chi)" 1 ) being higher in offshore surface 
waters, values of phytoplankton absorption spectra at the SeaWiFS wavebands were 
highly correlated to modeled values. The normalized scattering spectra! shapes and the 
mean spectrum were in agreement to observations in other coastal waters, while the 
backscattering ratios were on average lower in phytoplankton dominated surface waters 
(0.0101 ± 0.002) and higher in near-bottom waters (0.0191 ± 0.0045) with low Chi. 
Average percent differences in remote sensing reflectance R rs derived from modeled and 
in- water radiometric measurements were highest in the blue wavebands (52%) and at 
sampling stations with a more stratified water column. Estimates of Chi and CDOM 
absorption derived from SeaWiFS images generated using regional empirical algorithms 
were highly correlated to in situ data. 


OCIS codes: 010.4450 , 290.1350, 120.0280. 
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1. Introduction 


The Mississippi River (MR) system being one of the ton 10 rivers in length, 


freshwater discharge, and sediment delivery 1 , strongly influences primary productivity 

and carbon fluxes in the Northern Gulf of Mexico 2 . The large spatial and temporal 

variations in biogeochemical processes occurring in these river dominated coastal 

margins coupled with large seasonal variability of freshwater discharge, coastal 

circulation and effects of frontal passages or storms have contributed to general 

uncertainties in the knowledge of oceanic fluxes of organic and inorganic matter to the 

shelf and open ocean waters/ Satellite remote sensing with its capability for synoptic and 

repeated coverage of coastal waters has been found useful for monitoring many coastal 

processes. Ocean color remote sensing has been used to obtain estimates of 

phytoplankton biomass, colored dissolved organic matter (CDOM), water clarity, primary 

productivity and for tracing oceanic circulation features useful for carbon cycle studies. 4,5 

In recent years, occurrence of hypoxia or low dissolved oxygen concentration in the 

Louisiana coastal waters have been associated with high primary production and algal 

2 6 

blooms due to enhanced nutrient concentrations in waters discharged by the MR. ’ 
Reliable monitoring of these algal blooms will require improved knowledge of the optical 
characteristics of seawater constituents, and better optical models relating seawater 
constituents to the inherent optical properties (IOPs) of absorption and backscattering and 
the apparent optical properties (AOPs) such as remote sensing reflectance R rs . 

Empirical and semianalytic algorithms have been developed to obtain estimates of 
oceanic constituents and optical properties from ocean color data or remotely sensed 
reflectance data. These algorithms appear to perform well in oceanic Case 1 waters where 



phytoplankton and its derived material are the main constituents determining ocean 

color. 7 while their oerformance has been less than satisfactory’ in many optically complex 


Case 2 waters where other constituents (e.g., CDOM, detrital material or suspended 
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sediments) are present and may not covary with chlorophyll a concentrations (Chi). 5 
Empirical algorithms use statistical relationships to link measured reflectance ratios to 
Chi, 10 while semianalytic algorithms use a combination of equations derived using 
radiative transfer theory and empirical relationships (e.g., parameterization of the 
absorption coefficients in terms of Chi) to provide linkages between the AOPs such as 
reflectance, the IOPs (e.g., absorption and scattering) and seawater constituents. 1 1 
Development of semianalytic algorithms have allowed for a greater number of variables 
(e.g., Chi, CDOM absorption) to be determined from remotely sensed ocean color 
satellite data 5 ' 12,13 and have been mainly applicable in Case 1 waters where Chi is the 
main variable influencing optical properties. Using an approach similar to Case 1 waters, 
studies have also addressed the application of empirical and semianalytic solutions to 
coastal waters including river impacted coastal systems, 8,9, 14,1 3 where for example, 16 a 
linear relationship was shown to exist between the IOPs (absorption, backscattering) and 
the AOPs (diffuse attenuation coefficient, reflectance) in waters impacted by the seasonal 
MR discharge. Physical processes such as cold fronts and wind induced coastal 
circulation 17,18 can further influence biogeoehemical processes in the northern Gulf of 
Mexico. Frequent outbreaks of cold fronts in the winter months (generally from October 
to April) have been shown to affect water column physical properties (e.g., deepening of 
mixed layer) and air sea heat fluxes due to drier and colder air mass over the oceanic 
waters. 19,20 
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While these physical influences (frontal systems, wind induced coastal circulation) 
have been studied with respect to plume dynamics, physical effects and air-sea 
interaction, their influence on the bio-optical properties and remote sensing algorithms 
have not been addressed and is largely unknown. A field bio-optical study in waters 
impacted by the MR that coincided with the passage of a cold front and strong westward 
circulation currents allowed us to examine their effects on bio-optical properties, and 
remote sensing semianalytic and empirical algorithms. Using an approach similar to that 
for Case 1 waters, we examined the various parameterizations of the IOPs (absorption 
and scattering) in terms of Chi, and the correlations between radiometric derived R rs and 
those derived from field measurements of IOPs. Further, we evaluated seasonal 
variability (April and October 2000, and March 2002) in optical properties to obtain a 
better understanding of the role of various seawater constituents (e.g., phytoplankton, 
detrital or nonalgal particles, and CDOM) on the optical properties and ocean color 
algorithms in an optically complex river dominated coastal environment. 

2. Data and Methods 

A. Study site and sampling strategy 

The study site (Fig. 1) located on the Louisiana shelf in the Northern Gulf of Mexico is a 
region most directly influenced by the discharge from the MR. The MR enters the Gulf 
mainly through the birdsfoot delta through three large passes; the largest of which is 
Southwest (SW) Pass, followed by South Pass and East Pass. A seasonal river discharge 
and generally westward flowing alongshelf currents 17 results in a major region of 
freshwater influence along the coast of Louisiana and Texas. Frontal passages (between 
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months of October and April) associated with cold air outbreaks often pass through 


coastal Louisiana with characteristic time scales of 3 to 10 days. 


20 


Sampling was conducted aboard the R/V Pelican from 1 7 to 26 March 2002 (high 
discharge period) that coincided with the passage of a frontal system. Twelve stations 
were occupied during the cruise aligned along three North-South transects (Fig. 1) that 
extended from shallow (~ 8 meters) coastal waters to offshore stations on the continental 
shelf. To examine the influence of physical processes (frontal passage and coastal 
circulation), the data were partitioned into measurements taken before (17-21 March) and 
after (23-26 March) the frontal passage. All 12 stations were occupied before the frontal 
passage while all but stations 9 and 1 1 were occupied after the frontal passage. At each 
station, water column optical and hydrographic measurements were made using a bio- 
optical instrumentation package and a CTD instrument. Water samples were taken at 
discrete depths (usually three) using Niskin bottles and processed onboard the ship or 
later in the laboratory. Atmospheric conditions during the cruise were obtained from a 
NOAA National Data Buoy Center C-MAN station located near the SW Pass (Fig. 1). 


B. Hydrographic and bio-optical profile measurements 

A vertical profiling bio-optical package comprising of a CTD (SBE 49, Sea-bird 
Electronics) and a suite of optical instruments enabled simultaneous measurements of 
water column optical (absorption, scattering, basckscattering, chlorophyll fluorescence) 
and physical (temperature, salinity) variables. A data acquisition system (DFI4, WET 
Labs) was used to acquire and archive data collected from all instruments in the bio- 
optical package and subsequently transfer data to a personal computer. Time-stamped 
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data from these instruments were aligned to CTD depth and vertical profiles were binned 


at either 0.5 m (shallow water stations) or 1 .0 meter depth intervals. Profiles of 


chlorophyll fluorescence obtained using a flow through fluorometer (Wetstar. WET 
Labs) were converted to chlorophyll profiles by comparing High-Performance Liquid 
Chromatography (HPLC) Chi pigment estimates at discrete depths to fluorescence values 
at the same depths 21 and are based on the assumption of a constant quantum fluorescence 
yield. 22,23 Absorption and scattering coefficients at nine 10-nm spectral bands (412, 440, 
488, 510, 532, 555, 650, 676, and 715 iim) were determined using a 25-cm dual path- 
length absorption and attenuation meter (ac-9, WET Labs). Using optically clean water as 
a reference (obtained from daily calibration of the ac-9 during the cruise), profiles of 
absorption and attenuation were corrected for temperature, salinity and scattering"" 4 ’ 25 to 
give absorption a pg (particulate absorption a p plus a C dom) and attenuation coefficients c. 
Particulate scattering b p was determined as the difference of attenuation and absorption 
coefficients for all wavebands. Vertical profiles of absorption coefficients due to CDOM 
ciedom were obtained by attaching a 0.2 pm pleated maxicapsule filter (Gelman) to the 
intake tubing of the ac-9 during a separate cast at each station. 

Vertical profiles of backscattering coefficients at three wavebands (450, 532 and 650 
nm) were estimated using a factory calibrated backscattering meter (VSF-3, WET Labs) 
from measurements of optical scattering at three angles of 100, 125, and 150° at each 
waveband. The relationship of backscattering coefficient bb to the volume scattering 

♦ » "i 7 

function (VSF denoted by p(0)) through its integral in the backward direction"" is used to 
determine the particulate backscattering coefficient hb P at three wavebands from field 
measurements of particulate p p {p p (100°), p p (125°), and p p (150°)) respectively. The 
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measured j3 p (0) values were corrected for light loss due to absorption over the pathlength 
at each, angle and wavelength using the ae-9 data." 1 Further, a third-order polynomial was 
fitted through the three measured angles {2nf3 p (0)sin 6) plus a fourth angle 9 -tc, where 
f5 p (n)sin n~ 0. Particulate backscattering coefficient bb P at each of the three wavelengths 
were then determined by integrating the area under the polynomial. 

Remote sensing reflectance R rs (sr' 1 ) was determined from vertical profiles of 
downwelling irradiance Ed(X, z) (W nf nnf ) and upwelling radiance L U (X, z) (W nf nm" 

1 sr' 1 ) at seven 10-rnn wavebands (412, 443, 490, 510, 555, 66 5 and 683 nm) using the 
OCI and OCR-200 radiometers (Satlantic). Dark values at these wavebands were 
obtained from night casts of the bio-optical package and were found to be relatively small 
for the downwelling iiTadiance (~ 0.3 W m' nm" ) and upwelling radiance (~ 0. 1 W m‘ 
nm' 1 sr' 1 ) sensors, respectively. Spectral attenuation coefficients K ( j and Ki. w ? ere 
determined as the negative slope of In Ed or L u versus depth. An 1 1 -point moving linear 


regression of In Ed or L u versus depth was used to obtain vertical profiles of Kd and Ki. 
Using values of Kd and Ki calculated at the shallowest measured depths, the values of 
Ed(X,0~) and L u (X,0 ') were estimated just beneath the sea surface. These subsurface values 
were extrapolated to values just above the sea surface Ed(X,0 + ) and L U (X, (f ) using known 
values of Fresnel reflectance at the water-air interface (1.04, 0.54) 10,29 and were not 
adjusted for winds or sea state. R rs was then calculated as 


r„ = LufotfyE&cr). 


(i) 


In addition to the in-water derived values of R rs , remote sensing reflectance was also 
determined from above-water measurements. Measurements of sky radiance and reflected 
radiance from the water surface and from a plaque with known spectral directional 
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reflectance (10 % spectralon, Labsphere) were made from the ship deck using a fiber- 
optic based high resolution spectral radiometer (GER 1500, Geophysical and 
Environmental Research) following SeaWiFS protocol and data processed for 

29 30 

determination of above- water R rs . ’ 

SeaWiFS ocean color data were obtained from Goddard Distributed Active Archive 
Center (DAAC) and processed using NASA’s SeaDAS (version 4.8) software for 
estimating Chi using the OC4v4 algorithm and for R rs using the Standard SeaWiFS 


atmospheric correction^ 1 and the near infrared iterative method. 


32 


C. Discrete water sample analysis 

Water samples at discrete depths at each of the sampling stations were collected using a 
CTD rosette multi-bottle array (General Oceanic Model 1015) immediately prior to or 
following optical profiles and samples either processed at sea or taken to the laboratory 
for various types of analyses (e.g., pigment concentrations, particulate and CDOM 
absorption spectra). Filtered samples (one each for each depth at a station) were stored in 
liquid nitrogen and processed later in the laboratory for pigments and spectral absorption 
within six months of sample collection. Phytoplankton pigment concentrations were 
determined using the HPLC method'^ following standard recommended SeaWiFS 
protocol. 29 while total particulate a p and detrital or nonalgal absorption a ncip were 
determined using the quantitative filter pad technique (QFT). 34,3d Phytoplankton 
absorption a p h was then calculated as the difference between a p and a nap . Spectral 
absorption of CDOM was determined at sea by first filtering seawater samples through 
pre-rinsed 0.2 pm Nuclepore membrane filters followed by absorption measurements 
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36,37 


using a multiple pathlength aqueous capillar}'' waveguide based spectrophotometer 
within three hours of sample collection. Scdom was then calculated as the slope of the log 
of the CDOM spectral absorption between 370 and 500 am. Suspended particulate matter 
(SPM) concentration (mg f 1 ) was determined gravimetrically using standard 


procedures. 
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D. Bio-optical Modeling 

Models based on numerical simulations of underwater radiative transfer have been 
developed to relate AOPs such as subsurface irradiance reflectance R(A,0~) or the above 
water spectral remote sensing reflectance R rs to IOPs such as absorption a(X) and 
backscattering b(A) coefficients (mf ) of the medium through the approximation ’ 

R rs ~ 0. 54 (f/Q) x b b /(a + b b ) (2) 

where the value 0.54 accounts for the Fresnel reflectivities at the sea surface, /is a 
parameter that depends upon the solar zenith angle, the optical properties of seawater, 
and Q is a factor that relates subsurface upwelling irradiance to the nadir upwelling 
radiance (sr) (just beneath the surface) and are both wavelength dependent/ 1 The value of 
the ratio f/Q has been shown to vary less than /and Q individually for oceanic waters" 19 ' 40 
and were determined from field observations made in April of 2000. 16 These values were 
found to be 0.103, 0.109, 0.1 14, 0.1 10, 0.105, 0.101, 0.09, and 0.089 sr at the SeaWiFS 
wavebands of 412, 443, 490, 510, 532, 555, 650, and 67 6 nm, respectively. The 
absorption a(A) and backscattering b(A) coefficients (Eq. 2) are the sum of contributions 
by the various seawater particulate and dissolved constituents including pure water. 
Knowledge of the spectral IOP characteristics of the various seawater components 
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(phytoplankton, detrital or nonalgal particles, and CDOM), and their parameterization in 

terms of Chi have been used in bio -optical models to extract information (e.g., 

„ ^ 5II4I 

concentration of seawater constituents, lOPs) from ocean color data or R rs . 

The total absorption coefficient a (m' 1 ) is the additive sum of contributions by pure 
water a Wi 42 phytoplankton a p h , detrital or nonalgal particles a tiap , and CDOM absorption 
cicdom and can be written as 

a(A) = aJX) + a p h(A) + ttmap(A) “1“ Ucdom(A). (-f) 

In modeling total absorption in terms of Chi, phytoplankton absoiption a p h(A) has 
generally been expressed in terms of chlorophyll-specific absorption coefficient a* p h(A) 
(m mg Chf ) as 

ciph(A) = a p h( 440) a* p t,(A), (4) 

with a* p h(A) spectra having the a p h(440 )- normal i ze d spectral shape and a p f,(440) related 
to Chi through power law relationship given by 43 

a ph (440) - 0. 0403 [Chlf 668 . (5) 

However, variability in a* p h(/ i) has been shown to exist for different regions, 4 ,44 between 
species, 3 ”' and within species due to changes in growth irradiance and pigment packaging 
effects. 3 ' In a different approach, given a P h(440 % the a p h(A) spectrum was modeled using 
the empirical relationship of Lee et al. 45 expressed as 

ctph(V = {a 0 (A) + cii (A) l n[ a ph ( 440)] }a ph (440 ) , (6) 

where the values for a 0 (A) and aj(A) were empirically determined at 10 nm interval 
between 390 and 720 nm. The a p h(A) spectrum can thus be determined using Eq. (6) 
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given a value of a p h(440) that can be obtained from Chi. In this study we examine both 


approaches to determine a more appropriate model for the study region. 


With both CDOM and detrital spectral absorption exhibiting an exponential decrease 
in absorption with increasing wavelength, bio-optical models^' 41 have generally combined 
the contribution by the two seawater constituents to the absorption field using values for 
slope of the log-linearized absorption coefficient S in the range 0.014 to 0.020. These 
have been fitted to the model 


&cdom+ nap(A) ^cdom - nap(^-o ) C " " » (2) 

where a C d om +nap(Ao) the absorption at a reference wavelength (generally between 400 and 
440 nm) is then related to Chi. 40 In this study we examine the slopes S cdom and S„ ap of 
CDOM and nonalgal particle absorption coefficients separately using Eq. 7 to assess the 
appropriateness of using a single slope value in bio-optical models for the study, region. 

The total spectral scattering coefficient b(A) (m" ] ) can be separated into contributions 
by pure water b w 41 and total particulate b p as 


b(A) = b w (A) + b p (A). (8) 

With the backscattering coefficient being the term used in Eq. (2), it has, until recently, 
been a difficult parameter to measure directly. As such, particulate scattering b p has been 
related to particulate backscattering bb p through the backscattering ratio ( B = bb/b p ), and 
the backscattering coefficient expressed as 11 

btW= 1 /:b w (A) + Bb p (X).. (9) 

A constant value generally used for the backscattering ratio (B = 0.0183) was determined 
from field measurements conducted in both coastal and oceanic waters 48,49 Further, 
particulate scattering has been generally parameterized in terms of Chi as 30 
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( 10 ) 


b p (X) = 0.41 6[Chlf 766 (550/X). 


The spectral models for the lOPs and their parameterization in terms of Chi as 


described above have allowed the determination o+ seawater constituents and 10 Ps trom 
ocean color or R rs data mainly in Case 1 waters where phytoplankton and their derived 
products dominate the optical field. Other regional parameterizations of the IOPs have 
looked for similar dependence of absorption and scattering coefficients to Chi 46 Here we 
take a similar approach in examining the spectral characteristics of the IOPs and their 
relationship to Chi in efforts to develop bio-optical models for river dominated waters 
where Chi concentrations are often high and could play a major role in influencing 
optical properties. 


3. Results and Discussion 

A. Atmospheric conditions, circulation features, and surface water properties 
Prevailing winds in the MR delta region during March 2002 were easterly, southeasterly 
and have been shown to result in wind-driven currents that are primarily westward. 51 
These easterly winds are generally perturbed in winter (October to March) by frequent 
cold front passages that change the direction of the wind field from southeast to 
southwest, northwest, and northeast (Fig. 2). Relatively warm air temperatures (~ 22° C) 
with prevailing southeasterly winds (~ 6 m/s) during the prefrontal phase rapidly changed 
with the frontal passage to northerly winds (~15 m/s) resulting in a high pressure (1026 
mbar), cold air (~ 1 0° C) system occupying the study region. By 23 March, wind speeds 
were back to about 6 m/s although the cold air occupied the region for another 24 hours 
before rising again. 
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Based on the general easterly, southeasterly winds in the region, westerly circulation 


currents were expected 37 and also observed in the region from ADCP measurements (not 


shown). SeaWiFS Chi images (Fig. 3a, b, c, and d) documents the synoptic Chi pigment 


distribution before the frontal passage (F7, 19 March), during the frontal passage (22 
March), and after the frontal passage (24 March). These satellite images reveal the plume 
structure associated with higher Chi around the delta (Fig. 3a) and its intensification (Fig. 
3b) due to in situ growth. Before the frontal passage, higher concentrations of algal 
biomass associated with discharge from the SW, South and East Pass were concentrated 
around the delta (Figs. 3a, b) and extended south up to about 28.5° latitude into the Gulf 
of Mexico. High Chi associated with increasing river discharge from SW Pass (Fig. 2c) 
were observed in the satellite Chi images of 22 and 24 March (Fig. 3c, d) along with a 
separate S-shaped band of high Chi that extended from the nearshore (st. 8) to the outer 
shelf (approximately 28.5° south latitude) showing patterns of westward displacement. 
Previous studies 31 have documented a similar westward flow of plume waters turning 
anticyclonically northwestward towards the coast and have been attributed to the 
combined effects of volume discharge, winds, and the Coriolis force. 53 

Cross-shelf transects through stations 8, 7, 6, and 5 on 18 and 24 March 2002 (Fig. 4) 
provided spatial and temporal information on near-surface (~ 2m depth) bio-physical 
properties of temperature, salinity and chlorophyll fluorescence (obtained from the ship- 
based flowthrough system with a CTD and fluorometer) associated with physical 
processes in the region. Surface salinity on the cross-shelf transect on 1 8 March 2002 
(Fig. 4a) was relatively low (about 30.5 psu) except over a broad swath near st. 6 where it 
increased to about 32 psu. Phytoplankton biomass as indicated by the Chi fluorescence 
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signal varied from around 1 0 to 40 RFU (relative fluorescence units) and was lowest in 
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higher or lower with salinity. Salinity through the same transect on 24 March after the 


frontal passage was lower and more variable than on 1 8 March, indicating the coastal 
current transport of river plume waters across the transect (Fig. 4b). Temperature 
followed a similar pattern as that of salinity and generally decreased in lower salinity 
waters. The phytoplankton biomass was however highest near the offshore stations and 
lowest at the nearshore st. 8. A band (~ 5 km width) containing very high phytoplankton 
biomass (~ 120 RFU) was present between stations 5 and 6 and was almost six times 
greater than at the nearshore st. 8.(~20 RFU). The SeaWiFS Chi satellite imagery of 24 
March helps explain much of this short-term temporal and spatial variability in surface 
Chi fluorescence including the high Chi band between sts. 5 and 6 which appears to be a . 
part of the S-shaped band observed in Fig. 3c. Patterns in the synoptic SeaWiFS Chi ... s 
images however may not appear similar to the fluorescence transect due to strong 
westward coastal currents (speeds of about 1 m s' 1 were recorded by the ship-borne 
ADCP; not shown) and the 24-hours that was required to cover the transect between sts. 5 
and 8. 


B. Temporal changes in bio-optical properties 

We examine the short-term temporal surface properties (e.g., seawater constituents, 
physical, and optical) in relation to the frontal passage or the coastal currents at stations 
along the mid-transect (sts. 8, 7, 6 and 5), and two offshore sts. 4 and 12 (Table 1) using 
data acquired from the optical profiling package and discrete water samples. These 
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surface properties were an average of the top three meters, or in the case of discrete 

sampling within two meters of the surface. In some instances, due to variability in surface 


density stratification, differences were observed in the values obtained from discrete 


water sampling and those averaged from profiles. We observed that along the mid- 
transect and at the near-shore sts. 8, 7 and 6, salinities decreased, with the greatest change 
(3 1 .89 to 27.56 psu) occurring at st. 6. At the offshore sts. 5, 1 2 and 4, surface densities 
varied from a low of 20.8 to 24.1 kg m“ 3 mainly due to salinity increases associated with 
westward flowing currents moving oceanic waters over these stations. At the nearshore 
st. 8, SPM decreased while CDOM absorption increased slightly over the time period. An 
examination of absorption components at 443 nm indicated that percentage contribution 
by CDOM at this waveband increased from 37 to 64% over the six-day period. 

The backscattering ratio {B = b bp/bp) has been linked to the size of the nonalgal 
particles 54 (detritus and minerals) and to the refractive index with higher values of B 
associated with particles of higher refractive index (e.g., inorganic particles), and lower 
values associated with lower refractive index particles such as phytoplankton 
biomass. 2S ' :o ’ 36 We observed a general increase in B across the mid-transect (sts. 5, 6, 7, 
and 8), indicating a relative increase in the inorganic component, while at the two outer 
sts. 4 and 12, a decrease in B (at 532 nm) indicated a relative increase in phytoplankton 
biomass after the frontal passage. At st. 6, B increased from 0.0072 to 0.017, and was 
associated with a larger increase in bb P (0.0069 to 0.056 m" 1 ) than b p (0.962 to 3.350 m' 1 ). 
We observed that although Chi more than doubled, there was an order of magnitude 
increase in SPM concentrations indicating stronger correlations between phytoplankton 
biomass and b p and SPM and b hp , respectively. Similarly, while there was only a doubling 
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of a p h(443 ), we observed an order of magnitude increase in a nap (443) from 0.022 to 0.379 


m' 1 , indicating that a„ ap (443 ) is more directly correlated to SPM concentrations. The 
increase in B from 0.0072 to 0.017 also indicated the presence of higher concentrations or 


inorganic particles or smaller sized organic detrital particles transported by the coastal 
current after the frontal passage. High SPM concentrations at st. 6 also suggests that 
advective flows were able to transport large quantities of SPM over relatively large 
distances. A comparison of optical properties (absorption and scattering) at the three 
offshore stations 4, 5, and 12, shows the largest changes (decrease) occurred at st. 4 
which appeared to be occupied by more oceanic waters (salinity increased from 29.84 to 
34.35). Although. we observed a decrease in both Chi (4.57 to 0.99 mg m J ) and CDOM 
absorption (0.21 to 0.073 m' 1 ), a decrease in the backscattering ratio from 0.01 1 to 0.008 
(mean for 26 surface samples 0.01 1 ± 0.003) suggests a greater organic matter 
component in these oceanic waters. However smaller increases in salinity at stations 5 
and 12 were accompanied by smaller increases in the concentrations of seawater 
constituents and optical properties indicating river influences at these stations. The large 
variability in the IOPs, and concentrations of seawater constituents (Chi, CDOM, and 
SPM) could be explained in terms of a combination of factors such as river influences, 
coastal currents and biological influences such as growth. 

Water column physical and bio-optical properties before and after the frontal passage 
at sts. 6 (Fig. 5) and 2, 4, 5, 10 and 12 (Fig. 6) further demonstrates temporal changes in 
bio-optical properties over a six-day period. A large decrease in water column density 
(lower salinity and temperatures) following the front (Fig. 5) at st. 6 was mainly due to 
advection of lower salinity plume waters. As observed in Fig. 5a, the surface mixed layer 
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extended to about 3 meters with sigma-t then rapidly increasing from 22.2 to about 24.0 
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density with detrth to oceanic condition 


could be mainly attributed to salinity increasing to 36.2 at about 38 m. in contrast, water 
column density on 24 March at the same location was much lower (19.2 kg m' at 
surface) due to mainly advective flows of lower salinity waters in the region. Higher 
surface Chi generally decreased with depth (from about 2.1 and 7.2 mg m' in the surface 
to 0.2 and 1 .5 mg m° at 25 m depth) over the tw f o periods. On March 18, backscattering 
ratio ( B = b^/bp) was a low of 0.0072 in the surface mixed layer (Chi and SPM values of 
2.1 1 mg nf J and 0.54 mg l' 1 , respectively) (Fig. 5a) that increased with depth and 
decreasing Chi. Lower Chi (<0.2 mg m°) and higher SPM (2.98 mg l' 1 ) at 38 m resulted 
in a backscattering ratio ( B at 532 nm) of about 0.02 and could be mainly attributed to 
higher mineral or detrital component in the particulate pool. Similar inverse relationship 
between Chi and B observed in Case 2 waters were attributed to increasing refractive : 
index of particles in the water. 3 ’ In contrast, on 24 March the backscattering ratio was 
uniformly high throughout the water column (~ 0.017) in spite of higher Chi in surface 
waters (5.90 mg m°) and could be attributed to the lower salinity plume waters 
containing relatively larger concentrations of mineral and detrital particles throughout the 
water column (e.g., SPM concentrations of 5.67 and 4.35 mg l' 1 at surface and at 15 m 
depth). Although a plot of Chl/SPM ratio versus B (Fig. 7) shows a general increase in B 
with decreasing values of Chl/SPM ratios, in the case of station 6, B remained relatively 
constant in spite of Chl/SPM ratios decreasing from 1 .27x10' at surface to 0.64x 10' at 
1 6 m. One possible explanation for the constant B profile is that with st. 6 being 
approximately 62 km to the west of S W Pass, the relatively high concentrations of SPM 
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would mainly comprise of smaller size mineral and detrital particles being transported 

over longer distances in suspension. 

Higher concentrations of seawater constituents in the surface mixed layer were 
reflected in the higher values of lOPs (Figs. 5 and 6). Vertical profiles of CDOM 
absorption, particulate plus CDOM absorption and scattering were higher in surface 
waters and generally decreased with depth (Fig. 5c, d), with much higher values on 24 
March. CDOM absorption (shown at 440 nm) more than doubled while particulate plus 
CDOM absorption and scattering were almost factor of three larger throughout the water 
column. An examination of vertical profiles of salinity, Chi and the optical properties of 
absorption and scattering and B at sts. 2, 4, 5, 10, and 12 (Fig. 6a-j) obtained- before and 
itiafter the frontal passage also shows lower surface salinity values and higher Chi and 
^..absorption and scattering in near surface waters with the largest increase in absorption 
: : d^and scattering observed at station 10. Vertical profiles of B (Fig. 6e and j) shows a ; 

^-general increase in B with increasing depth. Average surface values for all the stations 
were 0.0101 ± 0.002 while the average bottom values were 0.0191 ± 0.0045 before the 
frontal passage that increased slightly to 0.0193 ± 0.0057 after the frontal passage. 

Changes in B due to the relative influence of various seawater constituents (e.g., 
phytoplankton, suspended sediments, and detrital material) and their size variations 
suggests the need for better characterization of the seawater constituents and their optical 
properties. Examining such relationships between optical properties and the seawater 
constituents have the potential to provide useful information on bulk particulate 
composition and as a tracer of water movement 28,55,56 in river-dominated coastal margins. 
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A seasonal comparison (April and October 2000, and March 2002) of CDOM 


absorption at 412 nm plotted as a function of salinity obtained from discrete water 


samples at different depths (Fig. 8a) shows an inverse relationship between CDOM 

absorption at 412 nm and salinity (r 2 = 0.95). In spite of the large seasonal river discharge 

3 -1 

during the three measurement periods (mean of 17,000, 5000 and 17,500 nr s 
respectively during April and October 2000 and March 2002) a conservative relationship 
is observed between CDOM absorption and salinity. High correlations are observed for 
CDOM absorption varying over two orders of magnitude in a study region that covered a 
larger area and seasonal river discharge in April and October of 2000 (see D’Sa and 
Miller 2003). A scatter plot of S cdofh (the slope parameter for the a cdom spectral decay) as a 
function of a cd om(412) (Fig. 8b) shows as in some previous results 57 S cdom decreasing with 
increasing a cdom (or decreasing salinity) for the March 2002 period while we observe a 
much greater scatter for the April and October 2000 cruises. Contributions of instrument 
variability to the scatter cannot be discounted since CDOM absorption measurements for 
the spring and fall 2000 cruises were made using an initial version of the multi-pathiength 
waveguide and a spectrophotometer, respectively. Average S cdom for April and October 
2000, and March 2002 varied over a relatively narrow range and were 0.015 1 ± 0.001 7, 
0.0162 ± 0.0018 nm' 1 , and 0.0161 ± 0.00054 nm' 1 , respectively. A closer examination of 
values of both CDOM absorption and S cdo m for the March 2002 measurements revealed a 
few stations with possible different sources of CDOM or composition such as humic and 
fulvic acids. 57 Station 9 had CDOM absorption values that were lower than expected and 
was associated with outflow from the Barataria Bay (measured surface salinities were 
lower than nearby stations 10 and 8). In comparison to values of 0.016 nm 1 at stations 2 


20 



and 3 observed before the frontal passage, S c d om values at these stations were 0.0149 and 
0.0148 nra' 1 or. 26 March and may have been related to a different source of river water 
associated with the increasing river discharge (Fig. 2c). However, the largest observed 
variation in S cc j 0 m (0.0141 nnf 1 ) occurred at station 6(18 March, 38 m depth, salinity 
36.14 psu) and appears to be associated with an anomalous water mass with slightly 
elevated Chi and scattering values (Fig. 6a, c). 


C. Characteristics of absorption coefficients 

IOPs and seawater constituents have been determined from inverse modeling of remote 
sensing reflectance using knowledge of the spectral characteristics of the surface water 
constituents and their relationships to Chi/ 1 5 We examined the spectral characteristics of 
particulate and dissolved surface seawater constituents (Fig. 9) and observed for example, 
a large range in the variability of spectral absorption values over a six-day period (pure 
water absorption shown for reference) for phytoplankton, CDOM and nonalgal 
absorption. For example in surface waters at station 6, phytoplankton and CDOM were 
the primary contributors to absorption before the frontal passage (Fig. 9a, a P hi ; Table 1), 
while CDOM and nonalgal particle absorption (a C( ] 0 m 2 , a nap 2) formed a significant 
component of the total absorption after the frontal passage. 

Phytoplankton absorption spectra obtained from surface samples at the sampling 
stations were highly variable and covered a large range at the blue (0.06 to 0.6 m' 1 ) and 
red (0.04 to 0.4 m' 1 ) absorption peaks (Fig. 9b, x-axis). Phytoplankton absorption spectra 
determined using the model of Lee et al. 43 with values of the constants ao(X) and aj(X) and 
Eq. 6 (at the SeaWiFS wavebands) were observed to be highly correlated to measured 
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phytoplankton spectra (Fig. 9b). However as observed in Fig. 9a (superimposed triangles 


and circles obtained using Eq. 6), the model fit at 555 nm 'waveband shows higher values 


than the actual a p h spectra, and as such by reducing the values of the two constants by 
0.0414 and 0.0061 (essentially using values of constants shifted up by 10 nm) we 
obtained a much improved fit (r 2 = 0.97; Fig. 9b) at this waveband for most a ph spectra of 
surface samples taken in the study region. 

Chlorophyll-specific phytoplankton absorption coefficient spectra a* p h of surface 
seawater samples (Fig. 9c) in the study area were observed to vary both spatially and 
temporally with variations in the blue and red absorption peaks of 0.02-0.1 and 0.02-0.09 
m 2 (mg Chi)' 1 , respectively. Highest values were observed at the offshore st. 5 for samples 
taken before and after the frontal passage and the iowest for stations directly influenced 
by the plume waters (sts. 2, 3, 6, 7, and 8). a* ph values for other stations were within 
those for sts. 2 and 5 . Factors such as pigment packaging, nutrients, taxonomic changes, 


and photoadaptation have been attributed to a* p h variability in oceanic waters 


35,43,58 


Many of these factors could have influenced a* p h variability in the study region which in 
conjunction with physical factors such as freshwater inputs, mixing and advection 
strongly influences primary productivity in the region. 2 Similar a* P h variability observed 
in other coastal waters 59 have important implications in the estimation of primary 
production from ocean color satellite data. 

In Fig 9a, we observed that although both detrital or nonalgal particles and CDOM 
absorption decreased exponentially with increasing wavelength, they do so with different 
slopes. Although most semianalytic models have combined the contributions by these 
two seawater constituents (Eq. 7), we observed that these slopes are different, being on 
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average about 0.01 1 for nonalgal particles (Fig. 9d) and 0.016 nm' 1 for CDOM (Fig. 9f). 


The a na p spectra of surface seawater samples followed the standard exponential decrease 
with increasing wavelength and the average slope compared well with those described in 
other coastal environments. 39,60 However there was a wide range in the magnitude of a nap 
absorption at 443 nm (0.019 to 0.892 m’ 1 ) being highest at the station near the SW Pass 
and lowest at the outermost oceanic stations. Although the average S nap slope was 0.01 1 
nm' 1 , it varied from 0.0085 to 0.0121 nm' 1 with the lowest values of 0.085 and 0.089 nm' 1 
observed at the nearshore stations 9 (surface waters) and 8 (bottom) off the Barataria Bay 
suggesting waters with different particle composition or properties. Absorption by 
nonalgal particles relative to total particulate absorption a p at 443 nm (Fig. 9e) varied 


from about 10% at the outermost. station 5 to about 85% for station 3 near the SW Pass 


and were in the range 10 to 30% for most near surface stations with chlorophyll w ' - 


concentrations between 2 to 6 mg. m°. The slope of a nap :a p versus Chi for stations 
sampled after the frontal passage appeared to be greater than those before the passage 
suggesting a greater nonalgal particle content in waters discharged by the Mississippi 


River after the frontal passage. Deviations were also observed in the CDOM spectral 
characteristics (Fig. 9f) where the CDOM absorption a cc i om (412) varied from 0.04 in 
oceanic waters to 1.2 m' 1 near the river mouth while the slope S C dom varied from 0.014 to 
0.0178 nm" 1 . These absorption properties suggest that near coastal waters with different 
freshwater sources potentially have optical properties different from surrounding water 
masses and the use of standard parameterizations may not be appropriate for these waters, 
In examining the absorption by phytoplankton, nonalgal particles at 443 nm and 
CDOM at 412 nm in terms of Chi (Fig. 10a, b and c), we observe that both a mp (443) and 
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cicdo,n(412) were highly correlated (r 2 > 0.9) to Chi for the March 2002 data while the 
April and October 2000 data showed greater scatter. Lower correlations between a P i,(443) 
and Chi (r 2 = 0.74) (Fig. 10a) for March 2002 data could be attributed to variability in 
surface Chi (standard deviation were ± 0.01 to 1.6 for Chi between 0.68 to 12.23 mg m* 3 ) 
resulting in some differences between discretely sampled values and Chi averaged over 
the surface 3 m. Comparison to limited data of April 2000 indicates similar relationship 
for the two periods. The non-linearity observed between suspended particulate matter 
(SPM) and Chi (Fig. lOd) for March 2002 were similar to those between a mp and Chi and 
could be attributed to the relatively large concentrations of nonalgal material (suspended 


sediments and detrital organic matter) at stations near the river. Strong currents observed 
A at the study site (>lm s' 1 ) may have carried large concentrations of SPM to distances 
4 such. as those observed at st. 6. However a nonlinear relationship is not observed for the 



il and October 2000 measurements with SPM that were low due to low river" •' ! 


4* discharge. We observe that the above parameterizations would not be applicable 
seasonally and also differed with those obtained Irom other studies, * suggesting a 
better need to better evaluate seasonal absorption parameterizations for retrieving in- 
water constituents and optical properties. 


D. Characteristics of scattering coefficients 

Relative dominance from contributions by different particulate field (e.g., phytoplankton, 
suspended sediments, detritus) influenced the spectral scattering coefficients b p which 
was highly variable (0.41 - 4.54 m" 1 at 440 nm; Fig 1 lb) and exhibited wavelength 
dependence that was rather small at lower absorption and increased at higher values (not 
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shown). Relatively different ratios of phytoplankton and nonalgal (suspended sediments 
and detrital) particles (Fig. 9e) before and after the frontal passage resulted in better 
correlations between phytoplankton and scattering before the frontal passage and between 
nonalgal particles and scattering after the frontal passage (not shown). River water 
discharge which increased dramatically after the frontal passage (Fig. 2c) appeared to 
contain higher concentration of nonalgal partic les that were better correlated to scattering 
(not shown). The scattering spectral dependence has been previously shown to vary as 2d 1 
with n generally having values from 0 to -2. 11,61 To account for the varying shape of the 
scattering spectrum, we normalized the b p (A) values with respect to the spectrally 


averaged scattering coefficient, <b p > that was computed between the 412 and 715 nm.' Q 
Troughs observed in the normalized scattering spectra are due chlorophyll a absorption 
effects on scattering corresponding to the absorption rnaxima in the blue (around 440 nm) 
and red (676 nm) spectral region; The normalized spectral shapes appeared similar to 
those in other coastal waters 09 with larger troughs at stations with high Chi and almost 


flat at locations with low Chi concentrations. The mean spectrum for all the data (Fig. 

1 1 a, bottom panel) compared well with the linear relationship between b p (X) and X for 
observations that included the northern Gulf of Mexico 49 and is flatter than the typical 
value of-1 used for the power function X n to describe the shape of the scattering spectra. 

Although a constant value (0.01 83) has been often used for the backscattering ratio 
5, 48,49 a nonlinear relationship previously reported for the study region 16 is being 
confirmed in this study (Fig. 9b). Water column profiles (Figs. 5a, b; 6e, j) and surface 
values of B at sampling stations (Table 1 ) confirm observations that B is not a constant in 
these waters, generally having lower values in phytoplankton dominated waters and 


25 



higher values in waters with higher concentrations of suspended sediments and indicative 


of the bulk refractive indices of the water. Increasing backscattering in relation to total 


scattering due to higher refractive indices in waters with higher suspended sediments or 
smaller sized nonalgal particles could be attributed to the observed non-linearity in the 
backscattering ratio. A comparison of the relationship between b p and Chi used in 
inverse modeling (e.g., Eq. 10) indicates that in these river dominated waters the 
observed relationship is higher than those developed for case 1 waters (Fig. lie). 30 As in 
other coastal waters* 9 a seasonal comparison of the data (spring and fall 2000) also 
indicates values falling above the model line that can be attributed mainly to increased 
contribution to scattering by the suspended sediments or detrital particles. 


E.lRemote sensing reflectance : radiometric versus model derived 
We compared the R rs spectra lor a set of representative stations derived from in-water 
radiometric measurements at the six visible SeaWiFS wavebands (412, 443, 490, 510, 
555, and 670 nm) with those modeled from surface IOP estimates (absorption, scattering, 
and backscattering ratio) applied to Eq. 2 (Fig. 12). We observed both high and variable 
concentrations of both SPM and CDOM in the study region. In such waters the bi- 
directional f/Q function besides being dependent on the Sun-angle, can be influenced by 

62 

both sediment and CDOM dominated waters or is simply unpredictable. However, we 
used the f/Q values obtained from a previous study in the region 16 on the assumption that 
the range of variability of the ratio//? is narrower than /or Q individually; 19,40 and the 
values could be typical for these coastal waters. 
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R rs spectra were found to be variable in shape and exhibited characteristics associated 


with oceanic waters where Chi was the main seawater constituent (st. 12; Fig. 12, third 


row) to more coastal waters (sts. 6, 10, 11,4, 7, and 8) where Chi, CDOM and suspended 

sediments all contributed to the reflectance field. The mean percent difference for all the 


Stations Shown in Fig. 12 (calculated as 100 X (R rS -modet ~ Rrs-msitu)/0'5(Rrs-m O del + Rrs-imito)) 
was highest in the blue at 412 nm (52 ± 28%) and lowest (-4 ± 18%) at the 555 nm 
waveband. However, at stations with low R rs values (e.g., st. 12 at 670 nm waveband) 


percent differences were large even though differences in absolute magnitude were small. 
Contributions by dissolved and nonalgal or detrital absorption significantly influenced R rs 
in the blue waveband at some stations, while backscattering by phytoplankton and 


suspended sediments influenced ocean color in the green and red parts of the spectrum at 
other stations. CDOM absorption at 412 nm which varied from a low of 0.103 (at st 12, 
sampled on 25 March) to a high of 0.417 m" 1 at st. 7 along with phytoplankton absorption 
(with Chi varying from 0:68 at st. 12 to 6.1 1 mg rn'lat st. 7) contributed more than 70% 


to total absorption in the blue wavebands and determined the R rs spectra with reflectance 


lower in the blue than the green wavebands. Relatively larger contribution (-30%) by 


nonalgal absorption (suspended sediments and detrital material) with SPM concentrations 


of 4.02 mg f 1 at st. 7 (bottom row) resulted in a large difference between the blue and 
green wavebands. The largest observed R rs value (0.00918 sr' 1 ) at the 555 nm waveband 
could be attributed to the high b hp (532) value of 0.0312 nf 1 measured at st. 7. The largest 


percent difference (107% in the blue and -40% in the red) between in- water and model 


derived R rs at this station may be attributed to the strong stratification where a surface 
low salinity (27.2 psu) mixed layer (0-4 m) with high Chi (6.1 1 mg m°) and a cc iom(412) of 
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0.416 rrf 1 overlay a higher salinity (30.3 psu) lower Chi (4.1 1 mg m’ 3 ) and CDOM 
subsurface water mass (Fig. 6. lower panel). Ignoring subsurface contributions of lOPs 
used in deriving the model estimates or the complexity of the underwater light field may 
have resulted in the large differences at this station. At st. 12 (third row) low Chi (0.68 
mg m' 3 ), CDOM (0.093 nf 1 ), SPM (0.98 mg l* 1 ) and backscattering b bp (532) (0.00536 nf 
] ) resulted in a shift of the peak in R rs from 555 nm to around 490 nm. Larger differences 
between the in- water derived and model derived estimates of R rs spectra in the blue could 
be attributed to the combined effects of higher light penetration in the blue and surface 
stratification, while stratification would have reduced effect in the red due to limited 


penetration of red light in the water. An examination of the vertical IOP profiles (e.gg 
Figs, 5, 6) indicated that the average surface (0 - 3 m) lOPs values used in deriving 'the 
Rrg § pectra using Eq. 2 may not al ways have been representative of their effects on the 
; radiance field. It appears that the highly stratified river dominated coastal waters with its 
no^Funiform vertical bio-optical properties may have influenced both the radiometric and 
IOP derived estimates of R rx . Smaller differences in R rs spectra observed at higher 


wavebands at most stations suggest that the model derived estimates are well represented 
by the IOPs in the surface 3 m, while in the blue contributions by CDOM and 
phytoplankton at greater depths influenced the radiance field. Also, since the 
measurements were made at different times during the day the upward radiance field 
could still be non-isotropic and dependent on Sun’s position in these coastal waters " and 
the f/Q values used may not have been the most appropriate. Other factors such as the in- 
water radiometric measurement uncertainties associated with extrapolation of subsurface 
measurements to above-water values, shelf-shading, and in the ease of modeled outputs. 
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to uncertainties in the scattering corrections for the ac-9 measured absorption and 


ittering coefficients 63 may have also contributed to the differences in R rx spectra. 

We also examined the SeaWiFS derived R rs spectra for three nearshore sts. 7 , 8 , and 


10 using the standard SeaWiFS atmospheric correction algorithm 31 that assumes 
negligible water leaving radiance in the near infrared pail of the spectrum and the 
iterative method 32 that corrects for the non-negligible water reflectance in the near 
infrared in waters with chlorophyll concentrations greater than ~2 mg m'\ We observed 
that at st. 7, the R rs spectra derived from SeaWiFS data using the two methods differed 
across the spectrum and had larger differences in the blue wavebands. The standard 


atmospheric correction algorithm compared well with the in-water derived R rs in' the blue 
i .-with both having low values due high absorption by C DOM, Chi and nonalgal or detrital 


..^particles, while the iterative algorithm overestimated the spectral values in the blue. In 
Tthe case of st. 8, a more uniform near-surface (0 - 6 m) high Chi (4.04 mg nW) and 
-4CDOM (0.312 m' 1 ) may have resulted in all the derived R rs spectra matching closely ' 
across the spectra. Both the SeaWiFS atmospheric correction algorithms correlated with 


the in-water and model derived reflectance suggesting the atmospheric correction 
algorithm to have worked well at this station. However at st. 10, R rs derived from 
SeaWiFS using both the atmospheric correction algorithms overestimated those derived 
using the model and in- water measurements (Fig. 12, second row). Although surface Chi 
(3.22 mg m' 3 ), CDOM absorption at 412 nm (0.29 m' 1 ) and SPM (3.76 mg 1 [ ) were 
typical of other stations, observed subsurface Chi (4.7 mg m°) peak (Fig. 6; top row) at 
about 8 m depth does not explain the differences between field and satellite estimates. 
The results of these comparisons suggest both the potential and complexity in developing 
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semianalytic and IOP inversion models and their application in ocean color remote 

sensing in these coastal waters. 

A comparison of R rs spectra derived from above-water measurements with those 

modeled using IOPs and the in- water radiometric derived spectra show the above-water 
derived R rs to have lower values across the spectra for all the stations. Better correlations 
between SeaWiFS, model and in-water derived R rs (Fig. 12) suggest the above-water 
measurements may be biased due to environmental factors such as Sun angle, cloud 
cover, wind speed, or viewing geometry 6 '. However, an examination of R rs spectra in the 
near-infrared (e.g., 670 11 m) indicates the above-water R rs to differ significantly from the 
other methods which appear to be in better agreement at this waveband. We note here 
that in deriving above-water R rs using SeaWiFS protocol 29, j0 sky radiance is removed, 
from the ocean radiance signal followed by the removal of any residual signal at 750 nm. 
Thercorrection at 750 nm is based on the assumption that the water- leaving radiance at 
750tnm is zero and that the reflected sky radiance has a white spectral shape f It appears 
that the residual correction also removed a component of the water leaving signal in the 
near- infrared resulting in values of R rs that were lower than those derived by other 
methods. Other factors such as environmental and viewing conditions could have 
contributed to uncertainties with above-water reflectance measurements. However, ease 
of measurements with the above-water radiometer has allowed for more routine use of 
this method in determinations of remote sensing reflectance for ocean color applications. 

F. Empirical ocean color algorithms for Chi and CDOM 
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Although we observed differences in magnitude between in-water and above-water 
derived R rs , a comparison of the ratios of R rs used in ocean color algorithms (e.g., 
R rs (490)/R rx (S55), or R rs (510)//?,.., (5 5 5)) for measurements made during three field cruises 
in the spring and fall of 2000 and March 2002 showed them to be highly correlated (Fig. 
13a, b). Based on these results and using a more extensive set of remote sensing spectra 
derived from above-water radiometric measurements made during the three cruises, 
regression analysis of log-transformed reflectance ratio R rs (490)/R rs (555) and Chi and the 
ratio R rs (5\0)/Rrs (555) and CDOM absorption (Fig. 13c, d) indicated high correlation 
between the R rs (k) ratios and Chi (r 2 >0.8) and CDOM absorption at 412 nm (r 2 >0.7). 
Power law equations for Chi and CDOM absorption derived from the regression analysis' 
were determined as . , . < • 

Chi =1 ,63(i?„(490) / i?„(555))' 2 ' 512 (11) 

• acddm(412) -.0.227(^, s (510) / R rs (55Sy 2m2 . . .(12) 

The R rs (490)/^n (55 5) reflectance band ratio algorithm for Chi and the R rs (5 1 0)/R rs 
(555) band ratio for CDOM absorption have been shown to perform relatively well for 
the study region. 16 Estimates of Chi and CDOM absorption obtained from SeaWiFS data 
using these empirical relationships show good agreement with a seasonal set of in situ 
data obtained in April and October of 2000 and March 2002 (Fig. 14). Although both 
seasonal and short-term effects were observed in the optical properties, the bio-optical 
relationships determined for these waters provided a sound basis for retrieval of these 
variables using the empirical algorithms. 
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4. Conclusion 


We examined bio-optical properties in a large river-dominated environment during the 


passage of a frontal system in March 2002 that also coincided with westward circulation 
flows and an increasingly large discharge from the Mississippi River. Optical properties 
were also compared to seasonal data obtained in the spring and fall of 2000 that 
corresponded to a year w'hen the discharge from the Mississippi River was very low. 
SeaWiFS Chi images obtained before, during and after the frontal passage in March 2002 
revealed high Chi around the delta, its intensification and subsequent westward 
displacement. Lower salinities observed during cross-shelf hydrographic measurements 
after the . frontal passage also indicated westward transport of plume waters: 'A ■ 
combination of factors that included transport, in situ growth (observed in satellite 
image), and a large increase in river discharge influenced the spatial and temporal 
variability in seawater constituents and the optical properties. A relative increase in u ; 
percentage of nonalgal absorption relative to algal absoiption after the frontal passage 
was reflected on an average increase in the backscattering ratio (0.0074 to 0.014) at 
stations directly influenced by plume waters while at the offshore stations it remained 
low (0.0085). Non linearities observed between a r , ap (443) and SPM and Chi could be 
attributed to the high concentrations of nonalgal particles (suspended sediments and 
organic detrital material) at stations more directly influenced by river water. Similar 
correlations were however weak or nonexistent for the April and October cruises. The 
spectral absorption and scattering characterization of seawater constituents such as 
phytoplankton, CDOM and nonalgal particles were found to be similar to those 
determined for other coastal waters/’ 7,49,59,60 We observed a non-linear relationship for the 
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backscattering ratio that can be attributed relative changes in particle composition. The 


observed short-term variability in algal to nonalga! absorption will relate to uncertainties 
in describing absorption in bio-opticai models 41 with consequences in obtaining reliable 

estimates of seawater constituents and optical properties from ocean color remote 
sensing. 

A comparison of R rs spectra derived from in- water radiometric measurements and 
those derived using the IOP model 51 indicated larger mean differences (52%) in the blue 
waveband where contributions by CDOM and phytoplankton absorption were most 
significant. There appeared to be a relative increase in the contribution by nonalgal 
absorption or backscattering 64 at some stations after the frontal passage associated with a 
large increase in river discharge. A combination of factors including deeper light 
penetration in the blue, water stratification, model derivation using only the near surface 
(0 $3m) values of IOPs and instrument measurement uncertainties contributed to the 
relatively larger differences in the blue waveband. In contrast differences at higher 


wavebands (green and near-infrared) were relatively small (4%) due to the greater 
absorption of light by pure water and its penetration to shorter depths. A comparison of 
model and in-water derived R rs spectra with those derived from SeaWiFS satellite data 
using the standard atmospheric correction algorithm and the iterative coastal algorithm 
indicated the best correlation to occur at a station with uniform near-surface Chi 
distribution while greater differences were observed at stations with stronger near surface 
stratification and non uniform vertical Chi distributions. This study has provided better 
insights into the variability of the ocean color reflectance measurements derived from 
both field and satellite data in relation to seawater constituents in an optically complex 



river-dominated system. Estimates of Chi and CDOM absorption derived from SeaWIFS 


ocean color data using the regional empirical algorithms were found to be highly 
correlated to in situ measurements acquired during both high and low MR flows. In spite 
of a conservative relation observed between CDOM absorption and salinity for 
measurements presented in this and other studies, 65 results of this study indicate both 
short-term and seasonal variability in the bio-optical properties suggesting need for more 
studies and caution in the use of semianalytic algorithms for this complex river 
dominated coastal environment. 
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FIGURE LEGENDS 


Fig. i. Study area showing the locations of the 12 stations (plus signs) occupied during a 
research cruise 17-26 March 2002 in the Northern Gulf of Mexico. Also shown are the 
three passes of the Mississippi River delta including the Southwest (SW) Pass through 
which most of the River water is discharged. Atmospheric data for this study was 
acquired from the NOAA C-MAN station located at the tip of the SW Pass. 

Fig. 2. Atmospheric and Mississippi River water discharge conditions during the study 
period in March 2002. (a) Wind speed and direction, and (b) pressure and air temperature 
Recorded at the NOAA National Data Buoy Center C-MAN station (BURL 1) (Fig. 1) 
vindicated the passage of a cold front system around 22 -March, 2002. (c) Mississippi River 
•water discharge recorded at the U.S. Army Corps of Engineers monitoring station at . 
ITarbert Landing, Mississippi. 


Fig. 3. (Color) Chlorophyll concentration derived from SeaWiFS data for the study area 
on (a) 17, (b) 19, (c) 22 and (d) 24 March 2002. The pluses indicate the location of field 
observations. Clouds are masked to black. 

Fig. 4. Cross-shelf transect through stations 5, 6, 7 and 8 (offshore to nearshore) 
conducted over a one day period on (a) 1 8 March and (b) 24 March 2002 showing the 
temporal variability in sea surface temperature (SST), salinity and Chi fluorescence 
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(relative fluorescence units, RFU) measured with a flow-through system onboard the R/V 

Pelican. 


Fig. 5. Water column profiles of physical and bio-optical properties at station 6 obtained 
on (a, c) 1 8 March and (b, d) 24 March 2002. (a, b) Comparison of vertical profiles of 
density, Chi, and ratio of particulate backscattering to scattering bt p :b p , (c, d), absorption 
(particulate plus CDOM, and CDOM) and particulate scattering. 

Fig. 6. Comparison of vertical profiles of salinity(a, f) and bio-optical properties such as 
Chi (b, g), absorption (particulate plus CDOM) (c, h), particulate scattering (d„ i), and 
ratio of particulate backscattering to scattering bi ?p :b p (e, j) for stations 2, 4, 5, 1 0, and 12 
. qbtained during periods before (top panel) and after (bottom panel) the frontal passage.- 
■ additional profile at station 12 obtained on 25 March is plotted to show representative 
daily variation in salinity and bio-optical properties. 

Fig. 7. The ratio of ChkSPM as a function of the backscattering ratio bb p (532):b p (532) for 
samples obtained at surface, mid-depth and bottom at the 1 2 stations before and after the 
frontal passage. 

Fig. 8. Results of CDOM absorption plotted as a function of salinity and the slopes of the 
log-linearized absorption spectra S ccPjm for seasonal data obtained in April and October of 
2000, and March of 2002. (a) Salinity versus CDOM absorption a cc / om (412) (m" 1 ) for 
discrete samples (surface, mid and bottom depths) obtained at the twelve stations during 
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March 2002. Data are also plotted for measurements obtained in April and October of 


2000 for eomoarison. Solid line is the linear regression fit to all data points, (b) S ca < om (nm 


) plotted as a function of a c dom(412) (m" 1 ) for April and October 2000, and March 2002. 


Fig. 9. Spectral absorption properties of particulate and dissolved seawater constituents in 
surface samples, (a) Comparison of phytoplankton, CDOM and nonalgal particle spectral 
absorption coefficients (a P h, a C( ] om , and a nap ) at station 6 on 1 8 and 24 March 2002 
showing the large variability in absorption properties. Pure water absorption is plotted for 
comparison, (b) measured versus modeled values of phytoplankton absorption at 
SeaWiFS wavebands, (c) phytoplankton specific absorption spectra a* P h(l)' y (d) log- 
linearized absorption spectra of nonalgal particles S nap (nm' 1 ), (e) ratio of absorption by 
non-aigal particles to total particulates at 443 rim versus chlorophyll concentrations; and 
-(f) log-linearized absorption spectra of CDOM before and after the frontal passage. 

Fig. 10. Comparison of (a) phytoplankton, (b) (nonalgal particles (at 443 nm), (c) 

CDOM absoiption (at 412 nm), and (d) suspended particulate matter (SPM) as a function 
of Chi (mg m°) for surface samples obtained in March 2002, and April and October 
2000. Linear regressions (solid line) were determined only for March 2002 samples. 


Fig. 1 1. (a) Normalized particulate scattering spectra as a function of wavelength (top 
panel), (bottom panel) average value of the normalized particulate scattering spectra of 
surface samples, and the average scattering spectra of Gould et al. 49 (b) Particulate 
backscattering bb P versus particulate scattering b p for surface waters at 532 nm, and (c) 
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particulate scattering b p { 555) versus Chi (mg m‘ 3 ) for surface samples showing regression 
line (solid) for March 2002. Also shown for comparison are the scatter plots for 
measurements obtained in April and October 2000. The dashed line shows the Morel 

model 50 for case I waters determined at 550 nm. 


Fig. 12. Remote sensing reflectance spectra R rs at SeaWiFS wavebands (412, 443, 490, 
510, 532, 555, and 670 nm) derived from in-water radiometric measurements (solid line, 
solid circle), and IOP model derived (dotted line, triangle) using Eq. 2. Also shown for 
comparison are the R rs spectra determined from above-water measurements using the 
GER- 1-50.0 fiber-optic based radiometer. Also plotted; are the SeaWiFS derived R rs 
.spectra using the Gordon and Wang J (diamonds), and the iterative J (squares) 
atmospheric -correction algorithms for , the nearshore stations 1.0, 7 and 8 obtained within 3: 
-hours of satellite pass.; Right panel shows the % difference between the IOP model and 
in-water radiometric derived estimates of R rs . 

Fig. 13. Comparison between the ratios of remote sensing reflectance R rs determined 
from measurements with above-water and in-water radiometers at same locations in April 
2000 and March 2002 for (a) ^(490)/7?„(555), and (b) tf„(510)/^(555). Results of (e) 
surface Chi (mg m' 3 ) plotted as a function of ratios of R rs at 490 and 555 nm wavebands 
(i?,,(490)/idX555)), and (d) surface CDOM absorption (nT 1 ) at 412 nm plotted as a 
function of R rs ratio R rs (5\0)/R rs (555)) for April and October 2000 and March 2002. 
Linear regression fit (solid line) was determined using all data points. 
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Fig. 14. A comparison between in situ and SeaWiFS derived Chi concentrations and 


"DOM absorption at 412 nm for data acquired in April and October 2000. and March 


2002. (a) Chi retrieved from ratios of SeaWiFS derived R rs using Eq. 11, and (b) ODOM 


absorption (at 412 nm) retrieved from ratios of SeaWiFS derived R rs using Eq. 12 .Solid 
line is the linear regression fit to all data points. 
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Table 1. Surface (0 to 3 m) physical and bio-optical properties for stations 8, 7, 6, and 5 
(mid-transect) and stations 4, and 12 (offshore) for the periods 17-21 and 23-26 

(brackets) March 2002 corresponding to the pre- and post frontal passage. Station 6 is 
marked in bold to highlight the station with the largest relative change in optical 
properties. 


Properties 
Temp. (°C) 


St 5 


Salinity (psu) 
Sigma-t (kg m' 3 ) 
Chi (mgnv) 
SPM (mg T 1 ) 
a ph (443) (m- 1 ) 
a nap (443) (m 1 ) 
a C dom(443) (m' 1 ) 
% a c dom(443) 

% a ph (443) 
b p (443) (m 1 ) 

b p (532) (m 1 )' 

b bp (532) (m- 1 ) 
b bp ( 5 3 2)/b p (5 32) 


8 

19784 (19.02)" 
30.26 (29.85) 
21.19(21.08) 
4.40(4.04) 
3.15 (1.56) 
0.268(0.114) 
0.067 (0.034) 
0.201 (0.276) 

37.1 (64.0) 
49.4 (26.4) 
2.076 (0.556) 
1.955 (0.534) 
0.0167 (0.006) 
0.0085 (0.011) 


St. 7 

18.83 (18.82) ' 
30.08 (27.18) 
21.31 (19.14) 
4.11 (6.11) 
1.38 (4.02) 
0.263 (0.202) 
0.056 (0.184) 
0.204 (0.213) 
38.5 (35.2) 
49.7 (33.3) 
3.049 (2.260) 
2.970 (2.094) 
0.0199 (0.031) 
0.0067 (0.014) 


St. 6 

20.46 (19.04) 
31.89 (27.56) 
22.29(19.95) 
2,1 1 (7.24) 
0.54 (5.67) 
0.103 (0.230) 
0.022 (0.379) 
0.155 (0.231) 

54.2 (27.2) 
36.7 (27.2) 
1.017(3.440) 
0.962 (3.350) 
0.0069 (0.056) 
0.0072(0.017) 


19.23 (19.60) 
30.28 (31.70) 
21.30(22.35) 
2.86 (3.48) 
1.80 (2.10) 
0.229 (0.340) 
0.066 (0.166) 
0.208 (0.277) 
40.8 (35.1) 
44.9(43.1) 
1.581 (2.336) 
1.498 (2.412) 
0.0125 (0.021) 
0.0083 (0.009) 


St 12 

21.30 ( 1 9.73) 

31.00 (33.56) 
21.37(23.49) 
2.10 (1.92) 
0.50 (1.86) 
0.078(0.099) 
0.020 (0.023) 
0.123 (0.093) 

54.1 (42.1) 

34,3 (44.7) 

0.684 (0.872) 
0.683 (0.849) 
0.0063 (0.007) 
0.009 (0.008) 


St. 4 

19' 94 (20.36) 

29.89(34.35) 

20.88(24.16) 

4.57 (0.99) 

1.47(1.19) 

0.238(0.120) 

0:065 (0.032) 

0.210(0.073) 

40.4(31.4) 

45:7(51.7) 

2.303 (0.902) 

2.220 (0.828) 

0.0255 (0.007) 

0.0114(0.008) 


St. 8 
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